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ABSTRACT 

We have identified an X-ray transient (hereafter M15 X-3) in the globular cluster M15 from an 
archival Chandra grating observation. M15 X-3 appears at an X-ray luminosity of 6 x 10^^ ergs s~^ 
with a spectrum consistent with an absorbed power law of photon index 1.51±0.14. The object is 
identifiable in archival Chandra HRC-I observations with an X-ray luminosity of 2 — 6 x 10^^ ergs s~^ 
and apparently soft colors, suggesting a neutron star low-mass X-ray binary in quiescence. We also 
observe it in outburst in a 2007 Chandra HRC-I observation, and in archival 1994-1995 ROSAT HRI 
observations. 

We identify a likely optical/UV counterpart with a (possibly transient) UV excess from archival 
HST data, which suggests a main sequence companion. We argue that M15 X-3's behavior is similar 
to that of the very faint X-ray transients which have been observed in the Galactic Center. We discuss 
several explanations for its very low X-ray luminosity, with the assumption that we have detected its 
companion. M15 X-3's uniquely low extinction and well-determined distance make it an excellent 
target for future studies. 

Subject headings: binaries : X-rays — stars: neutron — globular clusters: individual (NGC 7078) — 
accretion disks 



1. INTRODUCTION 

Low-luminosity X-ray sour ces in globular c l usters 
were identified with Einstein (|Hertz fc Grindlavl Il983l ) 
as a separate class of X-ray sources from bright 
(10^^<Lx<10^^ ergs/s) LMXBs containing neutron 
stars (NSs). HG83 suggested that these systems were pri- 
marily cataclysmic variables (CVs), with similar accre- 
tion rates to the LMXBs but potential energy wells 1000 
times shallower, but also including s ome NS LM XBs in 
quiescence (qLMXBs). However, Ve rbunt et al.l (^984) 
argued that the brightest low-Lx systems (Lx>10^^ 
ergs/s) were too bright to be CVs, and must be qLMXBs. 

Many low-Lx cluster sources have now been observed 
with the Chandra X-ray Observatory^ and the answer 
has been mixed; some of the brighter low-Lx systems 
have been ident ified with qLMXBs, some with CVs 
(e.g. in 47 Tuc, iGrindlav etalll200lf ). while others re- 
main unidentified. The brightest few \ow-Lx systems 
(5 X 10^^<Lx<10^^ ergs/s, generally transient) remain a 
puzzle for either explanation. These objects are more lu- 
minous than known CVs, and are brighter than qLMXBs 
emitting X-rays from heat stored in the crust between ac- 
cretion events. However, they are fainter than expected 
for X-ray binaries undergoing outbursts driven by the 
standard disk instability model (King 2000). Eight such 
systems may have been seen, with various X-ray tele- 
scopes (Einstein, ROSAT, Chandra)^ in several globu- 
lar clusters, but none have been studied in great detail 
(Heinke et al. 2008, in prep). Onl y one of them (Source B 
in NGC 6652. [Heinke et al"]|200ll ) shows strong evidence 
for being a quiescent low-mass X-ray binary, though its 
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accretion history is uncertain. 

Thirty very faint X-ray transients (VFXTs) have 
been identified, mostly in the Galactic Center region, 
having peak X-ray luminosities 10^^ < Lx < 10^^ 
ergs/s and quiescent at least a factor of 10 lower 
(iMuno et al. 2005b"; "Sakano et al.l l2005l : IWiinands et all 
2006; Degenaar & Wijnands 20081). Explanation of their 
very low outburst luminosities and low inferred time- 
averaged mass transfer rates (M^ 10~^^ ^o) 
quire new paths of binary evolution. These objects are 
now being well- studied with Chandra, XMM-Newton, 
and Swift through monitoring of the Galactic Center, 
where the majority have been seen. There is some over- 
lap between th ese VFXTs and the class of low- luminosity 
X-ray bursters (iCocchi et al.ll2001l : iCornelisse et al.ll2002l : 
iHands et al.ll2004l ), which indicates that at least some 
of these systems are low- mass X-ray binaries. Some 
low-luminosity X-ray bursters show extremely short X- 
ray outbursts (|Wiinands et al.l [20071 ). while at least one 
VFXT, XMMU J174716.1-28104 8, shows years-long ac - 
tivity at Lx =fewxlO^'^ ergs/s (|Del Santo et al.ll2007D . 
We are undertaking a program to search for and study 
VFXTs in globular clusters, which provide opportunities 
to study the optical, ultraviolet, and soft X-ray emission 
not available in observations of the Galactic Center. 

The globular cluster M15 (NGC 7078) is of particu- 
lar interest for X-ray and optical astronomy. Ml 5 is a 
massive cluster that shows signs of advanced core col- 
lapse, with a radial stellar density profil e showing a 
steep cusp to within the central a rcsecond (iLauer et al.l 
[T991I : iGuhathakurta et all [19961 : iSosin fc Kind I1997D . 
Ml 5 shows a central concentration of relatively un- 
der luminous mass, that may be composed of com- 
pact objects (such as heavy white dwarfs and neut ron 
stars, iDuh et all [19971 . l2003l : iBaumgardt et al.ll2Q03D o r 
an intermediate-mass black hole jGerssen et al ] [2002l) . 
Its high central density also makes it a likely lo- 
cation for stellar encounters, which should produce 
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Fig. 1.— Chandra images of M15, from 2000 ACIS-HETGS, 
2001 HRC-I, 2004 ACIS-HETGS, and 2007 HRC-I data. All data 
are smoothed with a gaussian of kern el 0.75 ^^ The dwarf nova 
(DN) identified by Hannikain en et al.l (|2005[ ) and M15 X-3 are 
marked. Readout streaks can be seen at position angles ~ 30 and 
~ 100 degrees from the saturated LMXBs in the HETGS data. 
The greyscale is chosen to emphasize M15 X-3 in each image. 

LMXBs, CVs, and millisecond radio pulsars. Stud- 
ies in the X-ray (requiring the use of Chandra)^ opti- 
cal, and radio have indee d identified two bright LMXBs 
(IWhite fc Angelinil [2QQlh . sev eral CVs (iCharles^et all 
I2QQ2I : iHannikainen et al 2005: "Di ebah et all I2007D and 
eight radio p ulsars (^nderson 199^. W e use a distance 
of 10.3 kpc (|van den^osch et al]T20Q6h and a column 
density of Nh = 4.6 x 10^^ cm-^ (|Janulislll99l ) in this 
paper. Here we report the detection of a very faint X-ray 
transient in the globular cluster Ml 5. 

2. CHANDRA DATA ANALYSIS 

We reduced all 6 Chandra observations of Ml 5 in the 
Chandra data archive^ (see Table [T]). All data were re- 
duced and analyzed using CIAO 3.4.1. For the ACIS 
exposures, we generated new bad pixel files, reprocessed 
the ACIS data removing pixel randomization, and fil- 
tered on grade and status. For all data, we searched for 
periods of high, flaring background, but found none.^ 

Several X-ray sources are apparent in the images 
( F i^ . [1]) , including the tw o bright LMXBs reported 
bv lWhite fc Angelinil (|200ll ), two faint X-ray sources (a 
d warf nova, DN, and plane tary nebula, PN) reported 
bv IHannikainen et al.l (|2005l ), and the transient reported 
here. After trying various detection strategies, we chose 
the following strategies, all on 2'x2' images centered 
on M15's core. To locate the bright LMXBs in the 
HETGS observations, we ran the CIAO detection pro- 

^ http:/ /cda. harvard.edu/chaser/mainEntry. do 

^ A smooth 60% rise in the background over the length of Obs_ID 
2413 is still small compared to the scattered X-ray background from 
the LMXBs. 
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Fig. 2. Combined 2001 Chandra HRC images of M15, in- 
cluding all good events (left) and only a soft subset (25<PI<120, 
right), binned by 2 and smoothed with a gaussian of kernel 0.75 ^^ 
Ellipses encircle sources detected by WAVDETECT in each image, 
which are labeled (PN=planetary nebula, DN=dwarf nova). 

gram WAVDETECT without an energy band restriction. 
We did this in order to ensure the bright LMXBs re- 
tained as many of their photons as possible, as the zeroth- 
order images of the two bright LXMBs show substantial 
pileup.^ To detect the other faint X-ray sources in Ml 5 
during the ACIS observations, we ran WAVDETECT in 
the 0.5-7 keV and 0.5-2 keV ranges. 

We reprojected the three 2001 HRC observations and 
merged them to increase sensitivity. We ran WAVDE- 
TECT on two images; the full HRC dataset, and a 
dataset filtered to include only softer photons by includ- 
ing only events i n pulse- invariant (PI) channels between 
25 and 120 (cf. ICameron etld] 120071 ). We binned the 
HRC data by a factor of two. We show both HRC im- 
ages of M15 in Fig. [1 

We chose WAVDETECT detection thresholds of le-6, 
a conservative limit which should produce much less than 
1 false source in each image. We identified no sources in 
the 2000 HETGS data besides the two bright LMXBs. In 
the combined, full-energy HR C image from 2001 we iden - 
tified two sources noted by IHannikainen et al] (|2005l ). 
identified as a dwarf nova (DN, their source C) and a 
planetary nebula (source D; PN K648), and determined 
accurate positions for the two LMXBs. In the 2004 data 
we identify a new transient (hereafter M15 X-3), as well 
as the DN. In the 2007 HRC-I data we again detect the 
new trans ient, along with the two bright LMXBs (as re- 
ported bv lHeinke etal]l2007[ ). 

Our WAVDETECT run on the filtered HRC data did 
not detect the PN or DN, but did detect a faint (8-count) 
source consistent with the position of M15 X-3 (Fig. [2]). 
No other, spurious, sources were detected in these de- 
tection runs, so we believe that we have indeed detected 
M15 X-3 in quiescence. Its detection in the filtered HRC 
data suggests that its quiescent spectrum is very soft, as 

^ Pileup is the recording of more than one photon as part of the 
same event during a single frametime, which leads t o distortion of 
the energy spectrum and loss of events (|Davisll200in . 
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is typically observed for NS LMXBs in quiescence. 

We aligned our astrometry with the best ra dio- 
determined position of AC 211 (|Kulkarni et al.lfl99Ql ) for 
the HRC frames, precessing from B1950 to J2000 using 
the NED precession calculator. We derive our most ac- 
curate position for M15 X-3 using the offset from AC 
211 in the 2007 HRC-I data, in Tabled Positions de- 
rived from the marginal HRC-I 2001 detection and 2004 
HETGS detection are consistent within 3a with this de- 
tection, but are subject to systematic errors (due to the 
marginal 2001 detection and effects of pileup on the de- 
rived position of AC 211 in 2004). Our listed positions 
for the other sources are derived from the offsets from 
AC 211 using the 2001 HRC-I data. 

We do not detect sources A and B from 
IHannikainen et all ([2005.) . though we do not utilize 
their psf-subtraction method. We n ote that the HRC 
X-ray source identified with a PN by IHannikainen et ahl 
(|2005f ) is undetected in the 2004 HETGS dataset, 
suggesting that it has a soft spectrum (though not as 
soft as M15 X-3 in quiescence) or is variable. We defer 
further analysis of the other X-ray sources in Ml 5 to 
future work. 
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Fig. 3. Top: Chandra 2004 HETGS X-ray spectrum (data and 
best-fit absorbed power-law model) of M15 X-3. Bottom: Residu- 
als to best fit. 



2.1. X-ray Spectral Analysis 

Transiently accreting NSs in quiescence are usually 
seen to have soft, blackbody-like X-ray spectra, often ac- 
companied by a harder X-ray component generally fit by 
a power-law of photon index 1-2 (Campana et al. 1998). 
The harder component is of unknown origin; an effect of 
continued accret ion, or a shock from a pulsar wind have 
been suggested (|Campana et al.|[l998f ). The blackbody- 
like component is generally understood as the radiation 
of heat from the NS surface. This heat is produced 
by deep crustal heating during accretion, and is radi- 
ated by the crust on a timescale of 10^ years, prod ucing 
a steady quiescent thermal NS luminosity (Brown et al.l 
119981 : ICampana et al.lfT998l : iHaensel fc Zdunikil990l ). At 
higher luminosities (above ^ 2 x 10^^ ergs/s), the spec- 
tra of faint accreting NSs are generally d ominated by the 
hard component (e.g. lJonker et al]l2004l ). The spectra of 
cataclysmic variables can often be represented by one or 
more hot plasm a components, e .g. the MEKAL models 
in XSPEC (Lied ahl et al.l [19951 ). or, with low statistics, 
a power-law of photon index 1-2. 

ACIS: We extract a spectrum, background and re- 
sponse files for M15 X-3 from the 2004 ACIS data, us- 
ing the CI AO specextract script. The spectrum can 
be well fit by an absorbed power-law of photon index 
r = 1.51 ± 0.14 at the cluster's Nh value (see Table gj 
Figure [3|). Adding a hydrogen- atmosphere neutron star 
atmo sphere model (the NSATMOS model, iHeinke et al.l 
l2006l ) does not improve the fit, but we can place con- 
straints on the luminosity in a NSATMOS component 
of Lx(0.5-10 keV)< 8.9 x 10^^ ergs/s. A hot thermal 
plasma model (MEKAL) can also describe the spectrum 
adequately, with a temperature of 14^6'' keV. 

HRC: HRC data is not amenable to standard spectral 
fitting, but hardness ratios can be computed and some 
spectral information extracted from comparison with the 
response matrix provided by the Chandra X-ray Center.^ 

^ littp://nedwww.ipac.caltecli.edu/forms/calculator.html 
^ See http://cxc.liarvard.edu/cal/Hrc/RMF/ 
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Fig. 4. — HRC-I hardness ratios for various power-law spectra 
(photon index F = 1 to 5, blue solid line) and blackbody spectra 
(kT=30 to 1000 eV, red dashed line) using the HRC-I response 
matrix. We show the effect of increasing the N^. columns from 
4.6 X 10^'-' to 10^^ or 10^^ cm~^ for the power-law with photon 
index 2 (magenta dotted line). M15 X-3's quiescent colors and 1-cr 
uncertainties are also plotted, suggesting a very soft spectrum. 

Only 13 counts are present within a 1 source region in 
the combined 2001 HRC data, of which ~3 are prob- 
ably background. Given these limitations, the ratios of 
counts in PI channels 1:100, 100:140, and 140:255 suggest 
a power-law photon index > 4 (Fig. [4j). This suggests 
a very so ft source, such as is often seen from quiescent 
LMXBs (|Jonker et al.ll2004[ ). although we caution that 
the statistics are very poor and the response matrix un- 
certain. 

We estimate the 0.5-10 keV X-ray luminosity of Ml 5 
X-3 in quiescence as 6 x 10^^ ergs/s (for a F = 2 power- 
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Fig. 5. — Background-subtracted lightcurve of M15 X-3 during 
the 2004 ACIS-S/HETGS observation. 

law) to 2 X 10^^ ergs/s (for a blackbody of fe T=llQeV). A 
hydrogen atmosphere model (NSATMOS, I Heinke et alJ 
1^06) with T=57 eV and 10 km radius provides a soft 
spectrum and comparable observed countrate, and a 
bolometric NS luminosity of 8 x 10^-^ ergs/s. Such soft 
X-ray sources would be hard to detect in zeroth-order 
HETGS images; in the 2000 HETGS observation, only 2 
or 1 counts are expected, given the powerlaw or black- 
body spectral models above. Therefore the nondetection 
of M15 X-3 in the 2000 data is not surprising. 

2.2. X-ray Variability 

To assess variability within the 2004 observation, we 
produced barycentered lightcurves of Ml 5 X-3 and ana- 
lyzed them using HEASARC's XRONOS software^ We 
show a lightcurve (binned at 3000 seconds) of M15 X-3 
in FiglSl A test on this lightcurve gives a probability 
of 0.02 that M15 X-3 is constant during the 2004 ob- 
servation. The power spectrum shows red noise below 
3 X 10~^ Hz, but no obvious periodic signals. 

3. ARCHIVAL X-RAY DATA ANALYSIS 

Several previous X-ray missions have had sufficient res- 
olution to separate M15 X-3 from the other 2 LMXBs in 
M15 (though not to separate the two bright LMXBs). 
We analyze archival data from ROSAT HRI & PSPC, 
ASCA SIS, EXOSAT HRI, and Einstein HRI observa- 
tions to look for evidence of M15 X-3. First we con- 
sider the constraints possible over the past 12 years from 
RXTE. 

3.1. RXTE All-Sky Monitor Analysis 

We used data from the RXTE All-Sky Monitor^^ to 
constrain M15 X-3's flux history. The known M15 
LMXBs (AC 211 and M15 X-2) together contribute 
roughly 4 x 10^^ ergs/s (2-10 keV), and their combined 
luminosity is seen to be fairly stable to within a factor of 
two, apart from short periods wh en Ml 5 is near the Sun 
and the data quality is poorer ( White fc A ngelini 2001). 
Thus, an outburst from Ml 5 X-3 would have to be of 
similar luminosity to be potentially discernible. 

^ http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xronos/xronos. html 
http://xte.mit.edu 




Fig. 6. — ROSAT images of M15. A weak source can be seen at 
a position consistent with M15 X-3, 22.6" south of the two bright 
LMXBs, in the two HRI images. The image scahng is adjusted for 
maximum visibihty of Ml 5 X-3. 

M15 X-3's observed Lx in 2004 might be the peak of 
an outburst, or the tail of an outburst with a higher peak 
Lx' Ml 5 is not conclusively detected by the ASM dur- 
ing April 2004, at a time when Ml 5 X-3 was active. A 
2-sigma limit to the flux of Ml 5 during the two weeks 
preceding this observation is 3 x 10^^ ergs/s in the 2-10 
keV band, assuming a Crab-like spectrum; this is consis- 
tent with the X-ray luminosity from Ml 5 over the RXTE 
era. A rough limit on the flux earlier in 2004, and over 
the RXTE era generally, is 6 x 10^^ ergs/s. This suggests 
that any outbursts from Ml 5 X-3 in the RXTE era did 
not significantly exceed 3 x 10^^ ergs/s, at least for more 
than a few weeks (but see below). 

Our 2007 HRC-I observation was taken 2 months after 
the highest fluxes ever recorded from Ml 5 by the RXTE 
ASM (50 mCrab,^^ or 1.5 x 10^'^ ergs/s). It is possible 
that these high fluxes were produced by another X-ray 
source besides the two previously known LMXBs, AC 
211 and M15 X-2. Our HRC-I observation finds M15 X- 
3 active, but not brighter than in the three other active 
episodes. Thus we do not clearly resolve whether M15 
X-3 or another source may have produced the high flux 
episode from M15. 

3.2. ROSAT Data Analysis 

Two lengthy ROSAT HRI observations of M15 exist 
in the HEASARC archive (PI: Grindlay), in addition to 
two shorter PSPC observations (we use only the on- axis 
one); see Table E The ROSAT XRT/HRI combination 
has a FWHM of 3", sufficient to resolve M15 X-3 from 
the two bright X-ray binaries in the cor e of Ml 5, and 
an ab solute positional accuracy of 5.3" (|Deutsch et al.l 
[19981 ). Although our detection algorithms do not iden- 
tify a source at this location, a weak source is visible 
to the reader in the two HRI images (see Fig. [6l left), 
with 165±21 and 35±9 counts above the local back- 
ground within 5" circles. Assuming an absorbed power- 
law spectrum of photon index 1.5 such as seen in 2004, 
we estimate Lx(1994, 0.5-10 keV)~ 7.4 x 10^^ ergs/s and 
Lx(1995, 0.5-10 keV)- 7.7 x 10^^ ergs/s. 

3.3. Other HEASARC Data Analysis 

For the remaining high-resolution imaging X-ray 
datasets (ASCA SIS, ROSAT PSPC, EXOSAT CMA, 

Noted on |htt p : / /xt e . mit . edu / xt e-anno . ht ml | 
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Fig. 7. — Postage stamp images of Ml 5 from six imaging ob- 
servations (ASCA, EXOSAT, and Einstein). The positions of both 
LMXBs and M15 X-3 are indicated with circles of radius 3". No 
boresight corrections have been apphed to these images. 

and Einstein HRI), we can provide only upper limits to 
M15 X-3's flux. The first step is to clarify whether the 
known LMXBs are those producing the detected source 
in prior observations. Figure [71 shows images from ASCA, 
EXOSAT, and Einstein images. In each case the blind 
pointing accuracy of the telescope is thought to be much 
better than 22 ''at 90% confidence, and in each case (see 
below) the detected source is consistent with the posi- 
tions of the known bright LMXBs, and not with Ml 5 
X-3. 

The background is in each case dominated by the point 
spread functions of the nearby bright LMXBs. For sim- 
ple, conservative estimates of the upper limits in each 
case, we use a simple procedure. We compare the count 
rates in 5" circles around the peak of the observed source 
(taken to be AC211 and M15 X-2) and the expected lo- 
cation of M15 X-3 (after boresighting to align the de- 
tected source with the point midway between the bright 
LMXBs). We do not subtract the expected counts from 
the wings of the LMXBs' point spread function, as the 
shape of those PSFs are not necessarily well known; this 
makes our limits very conservative. We compute the 0.5- 
10 keV Lx for the Ml 5 LMXB assuming a power law 
of photon index 1.51, and multiply by the ratio of the 
counts in the two circles. The limits derived are given in 
Tabled! 

ASCA: One ASCA observation is available. The po- 
sition of a bright ASCA SIS source can be centroided 
with in 12" after appropr iate corrections to the astrome- 
try (|Gotthelf et al.ll2000[ ). We find that the position of 
the M15 X-ray source, according to ASCA, is 21:29:58.88, 
+12:10:12.97, 12" NE of the positions of the known 
LMXBs, and thus consistent. 

EXOSAT CMA: EXOSAT CMA images from three 
epochs are available, each using the thin filter. The EX- 
OSAT CMA provides positions accurate to within 8" 
The position of the single detected source is within 8" of 
AC211 and M15 X-2 in 2 of the observations, and 12" to 
the NW in the third observation. 

^'^ http:/ /heasarc. gsfc.nasa.gov/docs/exosat/exosat_about. html 
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Fig. 8. — Timehne of sensitive X-ray detections and upper hmits 
for M15 X-3, from the ROSAT HRI and Chandra HETG and HRC- 
I imaging. The times of two useful HST observing epochs are also 
marked with squares. 

Einstein HRI: Einstein HRI images for two epochs 
are available . Einstein positions a re thought to be accu- 
rate to 3.2 " (|Grindlav et al.lll984D . and the detected M15 
source is consistent with the positions of either bright 
M15 LMXB. 

3.4. X-ray Conclusions 

The ROSAT and Chandra data require at least two 
substantial changes in flux; a decrease between 1995 and 
2000, and an increase between 2001 and 2004 (Fig. [8]). 
Consecutive X-ray observations tend to find the source 
in the same state, suggesting that typical on and off pe- 
riods for M15 X-3 are a few years in duration. Further 
high-resolution monitoring observations would be greatly 
useful in constraining M15 X-3's recurrence history. 

Our conclusion from the analysis of archival X-ray 
datasets on M15 is that if M15 X-3 has reached high 
Lxs of > 10^^ ergs/s, it has done so rarely. The twelve 
most sensitive X-ray observations over 29 years all con- 
strain M15 X-3 to Lx < 10^^ ergs/s. We consider three 
scenarios for M15 X-3's mass transfer rate (we assume 
a 10-km, 1.4 Mq NS accretor for each): a) the obser- 
vations in hand are representative, and it spends ~half 
of its time "on" at Lx = 6 x 10^^ ergs/s; for this case, 
M ^ 2 X 10"^^ Mo/year, b) M15 X-3 has undergone 
outbursts in the past 12 years, which we have happened 
to miss. The RXTE ASM coverage indicates no out- 
bursts brighter than ~ 3 x 10^^ ergs/s in the past 12 
years (apart from the 2007 brightening, ^3.ip . so we as- 
sume a maximum Lx = 3 x 10^^ ergs/s. The twelve most 
sensitive existing X-ray datasets constrain Ml 5 X-3 to 
Lx < 10^^ ergs/s; so we assume that X-3 is in outburst 
at Lx < 3 X 10^^ ergs/s only 1/12 of the time. This gives 
a mass transfer rate of M < 2 x 10~^^ M0/year. c) M15 
X-3 has undergone major outbursts, but they were more 
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than 12 years ago; in this case, we cannot constrain its 
mass transfer rate. 

4. OTHER WAVELENGTHS 

4.1. Optical photometry 

Ml 5 has been the target of numerous HST optical 
campaigns. M15 X-3's location 20'' from the center of 
M15 reduces crowding, but unfortunately places it out- 
side the field of view of most HST high-resolution instru- 
ments during pointings at the core. (The 0.1'' pixel size 
of the WFPC2 Wide Field chips does not adequately re- 
solve this crowded area.) Three epochs of imaging cover 
M15 X-3 with the Planetary Camera (PC) chip, or with 
ACS/WFC. Unfortunately, during the ACS observation 
epoch (2006), the orientation of the spacecraft placed 
bleeding columns from a nearby bright giant onto Ml 5 
X-3's position. Thus we remain with only two usable 
epochs, in only three bands: F336W (1994), and F555W 
and F439W (2002); see Fig. [9l and Table El 

Since the 2002 data were obtained with dithering, we 
constructed drizzled reconstructions of the F555W and 
F439W PC fields. In order to generate approximate 
color-magnitude diagrams (CMDs), we performed aper- 
ture photometry on the drizzled frames and on a com- 
bined frame of F336W exposures. We applied the stan- 
dard calibration for the STMAG system to the resulting 
magnitudes. We then applied an aperture correction to 
an aperture of 0.5". We approximately aligned the re- 
sulting color-magnitude diagrams (CMDs) to the John- 
son system by adding zeropoint shifts appropriate to a 
star at the main-sequence turnoff (MSTO); these shifts 
are neglible except for F439W, which was increased by 
0.60 mag. No color-dependent transformations were em- 
ployed. The resulting CMDs are shown in Figs. [10] and 
[TTl Given the use of aperture photometry without neigh- 
bor subtraction, we have not plotted the photometry for 
the region of radius 4" about the cluster center in or- 
der to produce cleaner fiduc ial sequences. Comp arison 
with the (V, B-V) CMD of lvan der Marel et all (120021 ) 
shows reasonable alignment of the magnitude scales in 
the vicinity of the horizontal branch and main sequence 
turnoff. 

We use the relative positions of M15 X-3 vs. AC211 
given in Table 3 to derive a 0.14" (2cr) error circle for 
M15 X-3's optical counterpart in the HST images. We 
have found a likely counterpart, within 0.05" (< la) 
of the expected position, that shows evidence of possi- 
bly transient UV excess associated with phases of higher 
mass transfer. The 2002 data show a clear detection in 
but only a marginal detection in B (Fig. [9]), mak- 
ing the B photometry of this object challenging. Since 
it was well below the automatic detection threshold in 
5, it had to be handled as a special case. The result- 
ing B — V color index puts the object well redward of the 
main sequence but with very large errors that include the 
main sequence (Fig. [TQ|) . As the 2002 HST data were 
taken only 8 months after the 2001 HRC observations 
that found Ml 5 X-3 in quiescence, it seems likely that it 
was also in quiescence in 2002. In contrast, comparison 
of the 1994 U and 2002 V data implies a clear UV ex- 
cess in 1994 (see Fig. [TT]), relative to the main sequence, 
consistent with a mass transfer event. This makes sense, 
considering that the 1994 HST image was taken only one 
month before a ROSAT observation that found Ml 5 X-3 





^ M 







Fig. 9. — Finding charts for candidate counterpart of Ml 5 X- 
3. The location of Ml 5 X-3 is bracketed by (red) dashes in each 
frame. Upper left, ACS/WFC F606W image from 2006 (affected by 
column bleeding); upper right, WFPC2 PC F555W image (2002); 
lower left, WFPC2 PC F336W image (1994); lower right, WFPC2 
PC F439W image (2002). 
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Fig. 10. Color-magnitude diagram for M15, in F439W ("B") 
and F555W ("V") filters, both from the same observation epoch 
in 2002. (The calibration is described in the text and includes 
a zeropoint shift of 0.6 mag in F439W to bring it closer to B.) 
The proposed counterpart to M15 X-3 is indicated by the red dot 
with the large error bars to the red of the main sequence (though 
consistent with it). 

to be active (although M15 X-3 could have varied during 
that month). 
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Fig. 11. Color-magnitude diagram for M15, in F336W ("[/") 
and F555W ("y") filters. (The calibration is described in the text.) 
The proposed counterpart to M15 X-3 is indicated by the red dot 
with error bars to the blue of the main sequence. Note that the 
F555W and F336W frames were taken in 2002 and 1994 respec- 
tively. 

4.2. Radio images 

The position of M15 X-3 lies in the field of view of 
four VLA obse rvations reported in t he literature. A 1.4 
GHz image by iKulkarni et al.l (|l99Ql ) shows no emission 
from M15 X-3's location, for an upper limit of 0.3 mJy 
(4cr) on October 11-12 1988. A deeper observation at 
1.4 G Hz between April 9-14, 1990, by iJohnston efall 
(|l99ir ) allows an upper limit of ^1 30 /iJy (3a). An 8. 4 
GHz map from Feb. 11, 1991, bv iKnapp etld] ((1996!), 
finds a marginal 3cr detection (~70 /iJy) at a position 
consistent with Ml 5 X-3 (within 2"; the VLA was in 
hybrid C/D configuration, with a larger beam). A fourth 
VL A observati o n on Oct. 13, 2004 at 8.6 GHz, reported 
bv iBash et all (|2008f ). allows an upper limit of 25 /iJy 
(4a). 

5. DISCUSSION 
5.1. Nature of the accretor 

M15 X-3's position near the center of M15 strongly in- 
dicates that it is a cluster member. M15 X-3's X-ray lu- 
minosity measurements fall within a range that is reached 
by both quiescent LMXBs and cataclysmic variables. We 
provide three arguments for why Ml 5 X-3 is more likely 
to be an accreting neutron star: 

Strong X-ray variability: M15 X-3's X-ray fiux de- 
clined by two orders of magnitude between 1995 and 
2000, and increased again to its prior level between 2001 
and 2004. Such strong X-ray variability is a hallmark of 
LMXBs, but is rar e among most cataclysmic variables. 
iBaskill et al.l (|2005l ) observe that nonmagnetic CVs vary 
by one order of magnitude or less between observations, 
and do not reach maximum X-ray luminosities above 
3 X 10^^ ergs/s. The X-ray luminosities of polars do vary 
by two orders of magnitude on appropriate timescales. 



but t hey do not reach above 10^^ ergs/s (|Ramsav et al.l 
l2004f ). Intermediate polars do attain Lx > 3 x 10^^ 
ergs/s, but they rarely sh ow substantial long-term vari- 
ability (|Muno et al.ll2004D . 

X-ray to optical flux ratio: The lack of simultaneity 
complicates this comparison. However, the 1994 ROSAT 
HRI observation which found Ml 5 X-3 in its high state 
was only one month away from the HST WFPC2 F336W 
observations discussed in §4.11 Considering that M15 
X-3 seems to remain in each state for at least a year, 
we assume that these values can be compared. M15 X- 
3's V magnitude in 1994 must have been > 21 (oth- 
erwise it would not be blue during outburst, given the 
U constraint). Using \og{Fopt) = —OAV — 3.96, we 
find Fx (0.5-2.5)/ Fopf=0.p , which is higher than for most 
CVs (jVerbunt et al.l1l997h . M15 X-3's relatively high X- 
ray to optical fiux ratio suggests a neutron star accretor. 

Variation in X-rays vs. optical: To within a large un- 
certainty in the B magnitude. Ml 5 X-3 is consistent with 
being on or near the main sequence in the (F, B — V) 
CMD, and thus is unlikely to be much fainter in qui- 
escence. Our outburst U detection thus limits the to- 
tal V increase in outburst to one magnitude, while the 
X-ray fiux varies by two orders of magnitude. Most 
CVs show larger variations in optical than X-ray (e.g. 
Garnavich & Szkody 1988). 

These arguments suggest that Ml 5 X-3 contains an 
accreting neutron star, although confirmation (through, 
e.g.. X-ray spectroscopy in quiescence) is desirable. If 
M15 X-3 contains a neutron star, the observed outbursts 
are among the faintest seen from neutron star accretors, 
identifying M15 X-3's behavior as that of a very faint 
X-ray transient. 

5.2. Very faint X-ray transients 

Very faint X-ray transients (VFXTs) are identified 
as X-ray transients with peak X-ray luminosities of 
Lx = 10^^-^"^"^ to 10^^ ergs/s a nd quiescent lumi- 



nosities at least 10 times fainter ( Muno et al.l l2005bl : 
ISakano et aP l2005l : I Wiinands e^dH l2006[ ) . Thirtv sys- 
tems showing such behavior have been identified, with 
3-4 of them being active at any moment in the Galac- 
tic Center. Some rise up to Lx ^ 10^^ ergs/s, with a 
low (< 10%) duty cycle, while others, such as XMMU 
J174716.1-281048 and CXOGC J174535.5-290124, main- 
tain Lx ~ iq34 ergs/s with rel atively high (^ 50%) 
duty c ycles (Muno et al. " 2005bl : iDel Santo etall 120071 : 
iDegenaar fc Wiinands 20081). M15 X-3 appears to show 
similar behavior to the latter group, being observed to 
be "on" at Lx ~ 6 x 10^^ ergs/s in half of the X-ray 
observations with sufficient sensitivity to detect such a 
state. 

Typical, better-understood X-ray transients have peak 
X-ray luminosities of Lx = 10^^ to 10^^ ergs/s, and time- 
averaged mass transfer rates of M>^10~^^ M0/year. The 
VFXTs, by contrast, show combinations of duty cycles 
and peak L^s which indicate time- aver aged mass trans- 
fer rates one to two orders of magnitude lower than 10~^^ 
Mrr^/yea r, difficult to understan d in a binary evolution 
context (|King Wiinandsll2006l ). We note that VFXTs 
may undergo normal outbursts at other (so far unob- 
served) times, and thus might not have particularly low 
time-averaged mass transfer rates. In either case, their 
unusually low observed accretion luminosities are diffi- 
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cult to understand within the standard disk instabihty 
model, and thus warrant study. 

Several possible explanations for VFXT behavior can 
be immediately ruled out for M15 X-3 (though they may 
be relevant for other VFXTs). Accretion fro m the wind 
of a massive (e.g. B) star or a subgiant (|Pfahl et al.l 
[2002) can be ruled out for an optically faint star in a 
glob ular cluster. An interme diate-mass black hole accre- 
tor (|King fc Wiinandsl[2QQ6h can also be ruled out, due 
to M15 X-3's offset from M15's core, and the expected 
rapid settling of massive objects into the core. 

If our observations are representative, we can infer 
a time-averaged mass transfer rate of M ~ 2 x 10~^^ 
M0/year ( §3.4p . This is roughly the minimum mass- 
transfer rate achievable through binary evolution of a 
main-sequence star with a neutron star co mpanion af- 
ter a Hubble time (|King fc Wijnandsll2QQ6l ). suggesting 
an extremely low-mass brown dwarf companion might 
be a possible solution. However, our HST photome- 
try suggests (with large uncertainties) that the compan- 
ion to M15 X-3 is a main-sequence star with a mass o f 
~ 0.65 Mq (using an isochrone from iBrilev et al.ll2QQj ). 
Confirmation of a main-sequence companion (through 
deeper high-resolution HST imaging) would clearly rule 
out an extremely low- mass substellar companion. It 
would also severely constrain possible models of accre- 
tion from the wind of a low- mass main sequence star, as 
IWillems fc Kdbl ([2003h found a maximum time- averaged 
bolometric accretion luminosity of 10^^ ergs/s from stars 
in this mass range in the pre-LMXB stage. 

One plausible alternative for M15 X-3 is an extremely 
unfavorable i nclination angle, suc h as the X-ray transient 
identified bv iMuno et al.l (|2QQ5al ). allowing an isotropic 
X-ray luminosity higher by two orders of magnitude. At 
10 kpc, an outburst reaching 10^^ ergs/s might be pre- 
dicted to reach radio fluxes of 20 /iJy from a neutron 
star LMXB, or 1.0 mJy from a black hole LMXB in the 
hard state (using the correlation from iMigliari fc Fenderl 
[2OO6). The stringent Oct. 2004 radio upper limit of 25 
/iJy, 6 months after an X-ray observation showing ac- 
tivity, suggests that if M15 X-3 is highly inclined, it is 
more likely a neutron star LMXB. The possible 70 /iJy 
radio detection in 1991, if real, suggests either a major 
1991 outburst (unfortunately not covered with sensitive 
imaging), or that a large fraction of M15 X-3's accre- 
tion power at low accretio n rates is going in to the for- 
mation of jets (as found bv lGallo et al.|[2QQ6l . for black 
hole systems). A high-inclination explanation for the low 
Lx might predict detectable eclipses (the 2004 lightcurve 
(FiglU excludes eclipses longer than 1000 s), and that 
the neutron star should not be directly visible in qui- 
escence as a blackbody-like X-ray source (testable with 
ACIS spectroscopy in quiescence). 

A second alternative is the removal of most mass 
transferred from the companion out of the system, per- 
haps by the pressure of a fast-spinning NS magneti c 
field, the "propeller effect" jm arionov fc Sunvaevl[l 975h. 
Detailed simulations of the interaction of an accretion 
disk with such a spinning magnetic field find that some 
material does leak dow n onto the neutron star poles 
(|Romanova et al.l l2005l ). A moderately effective pro- 
peller effect may explain the particularly low mass trans- 
fer rates inferred for VFXTs. In this case, the divid- 



ing line between VFXTs and "normal" LMXBs may 
be ambiguous, dependent on spin rates, magnetic field 
strength and/or geometry. VFXTs may slowly accumu- 
late sufficiently massive disks to produce normal out- 
bursts. Some objects do show VFXT-like outbursts as 
weh as brighter outbursts, e.g. GRS 1741.9-2853 and 
XMM J174457-2850.3 (Muno e£aD l2003l : ISakano et all 
I2OO5I : [Pegenaar fc Wijnandsn 20081 ). supporting this pos- 
sibility. 

A third alternative is variability in the mass transfer 
rates from the companion, on timescales greater than 
tens of years, which is thought to oc cur in cataclysmic 
variables (e.g. iHameurv et 311119891 ). (|King fc Wijnandsl 
120061 ) dismiss mass transfer cycles as an explanation for 
VFXTs as a class for two reasons: a lack of systems 
much brighter than the mean mass transfer rate to bal- 
ance those which are fainter, and a lack of a mechanism 
that can deliver such strong variability. However, the 
known bright neutron star LMXBs tend to be brighter 
than evolutionary calculations predi ct at their orbi tal pe- 
riods (Podsiad lowski et al. 2002; P fahl et al.ll2003r ). One 
very plausible driver of such variability is an irradiation- 
driven instability, wherein the X-ray irradiation of the 
secondary causes slow expansion, le ading to mass trans- 
fer cy cles on ^10^ year tim escales (Podsiadlowski" 199!; 
.Hameurv et al.lll993l : iRitter et al.,.2000. : .Biining fc Ritteri 
200I . 

6. CONCLUSIONS 

We have identified a transient X-ray source in the glob- 
ular cluster M15, denoted M15 X-3. This source has been 
detected at = 4 - 8 x 10^^ ergs/s in 4 ROSAT and 
Chandra observations, and Lx ~ 2 — 6 x 10^^ ergs/s 
in 3 Chandra HRC-I observations, with an apparently 
softer spectrum. Archival HST observations reveal a star 
with an apparent ultraviolet excess at the X-ray position. 
Within broad uncertainty, due to a very weak detection 
in 5, the location of the proposed counterpart in the {V ^ 
B — V) CMD is consistent with it being a main-sequence 
star. The X-ray and optical characteristics of this source 
suggest an accreting neutron star nature, with the X- 
ray flux history resembling those of "very faint X-ray 
transients" seen in the Galactic Center region. A 1991 
VLA observation may have detected Ml 5 X-3 at ~ 70 
/iJy. If M15 X-3 consists of a neutron star accreting from 
a main-sequence star, as the evidence suggests, the ob- 
served accretion rate onto the NS must be rather lower 
than the expected mass-transfer rate for a Roche-lobe 
filling system. 

The location of this object is a mixed blessing. Be- 
ing located in M15, 20" from two luminous and persis- 
tent X-ray binaries, hides this object from all X-ray in- 
struments but those with the highest angular resolution 
(only Chandra for the foreseeable future). Being in such 
a dense globular cluster also prevents optical followup 
with any instrument with inferior angular resolution to 
HST. However, M15's intrinsic interest has provided us 
with numerous archival observations that allowed the dis- 
covery of Ml 5 X-3. Most importantly, the low extinction 
in this direction allows future study of this object in qui- 
escence, which is excluded for all other known very faint 
X-ray transients both in the optical/UV and soft X-rays. 
M15 X-3 may be a key object for understanding the na- 
ture of very faint X-ray transients. 
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TABLE 1 

Summary of X-ray Observations of M15 X-3 





WDb. UdAiG 


rL/xposure 


Instrument 




4879-4882 


1979 Nov 19 


5.07 


Einstein-HRI 


<3.0e34 


656 


1978 Nov 22 


0.26 


Einstein-HRI 


<le35 


657 


1979 May 17 


1.60 


Einstein-HRI 


<3e34 


28831 


1984 Jun 30 


23.73 


EXOSAT-CMA 


<3e35 


37046 


1984 Oct 22 


36.15 


EXOSAT-CMA 


<9e34 


63182 


1985 Oct 20 


20.27 


EXOSAT-CMA 


<1.6e35 


400081N00 


1992 May 16 


8.79 


ROSAT-PSPC 


<3.4e34 


400611N00 


1994 Nov 27 


52.47 


ROSAT-HRI 


7.4± 1 X 1033 


400611A01 


1995 May 11 


27.62 


ROSAT-HRI 


7.7±2e33 


43042000 


1995 May 16 


37.38 


ASCA-SIS 


<2.7e36 


675 


2000 Aug 24 


20.03 


CXO-HETG^ 


<1.8e32 


1903 


2001 Jul 13 


9.10 


CXO-HRC-I 


2 X 1031 


2412 


2001 Aug 03 


8.90 


CXO-HRC-I 


2 X 103i 


2413 


2001 Aug 22 


10.87 


CXO-HRC-I 


2 X 1031 


4572 


2004 Apr 17 


60.14 


CXO-HETGS 


5.9 X 1033 


9584 


2007 Sep 05 


2.14 


CXO-HRC-I 


4.5 X 1033 



Note. — Times in kiloseconds. Unabsorbed Lx (0.5-10 keV) observed, or upper limits 
reached (see text), in each observation. We use PIMMS and a powerlaw of photon index 1.5 
absorbed by Nfj = 4.6 x 10"^*^ cm~^ to estimate detections or upper limits for observations 
where Lx may be > 1033 ergs/s, or PIMMS and a blackbody of 0.135 keV with the same 
for observations where Lx is constrained to < 1033 ergs/s. 



TABLE 2 
Summary of HST Observations 





Prog. ID 


Epoch 


Instrument Exposures 


Magnitude 






5742 


1994 Oct 26 


WFPC2 3x600 F336W 


21.5±0.2 






9039 


2002 Apr 5 


WFPC2 12x16 F555W 


22.0±0.2 






9039 


2002 Apr 5 


WFPC2 4x40 F439W 


23.7±0.8 






Note. — Maj 


gnitudes refer to proposed M15 X-3 counterpart, with calibration as discussed 






in text. 
















TABLE 3 










Chandra Source List 






Name 


RA 


Dec 


Cts: 2000 Cts: 2001 


Cts: 2004 


Cts: 2007 


AC211 


21:29:58.310 


12:10:02.85 


80769(299) 




4735(70) 


M15 X-2 


21:29:58.124(1) 


12:10:02.37(1) 


191255(454) 




13960(120) 


DN, src C 


21:29:57.323(6) 


12:10:43.9(1) 


[5] 21(5) 


79(10) 


< 1 


PN, K648 


21:29:59.40(1) 


12:10:27.6(1) 


< 2 9(3) 


[5] 


< 1 


M15 X-3 


21:29:58.155(3) 


12:09:40.04(6) 


< 2 8(3) 


480(23) 


31(6) 



Note. — Positions and counts from WAVDETECT, in 0.5-7 keV band for ACIS data, full band for HRC data. Nondetections 
indicated with upper limits, subjective detections in brackets. Errors are in parentheses, on the last quoted digit, and include 
only random (not possible systematic) errors. The exposure times for the observations are: 2000 ACIS-S/HETGS, 20.03 ks; 2001 
HRC-I, 28.87 ks; 2004 ACIS-S/HETGS, 60.14 ks; 2007 HRC-I, 2.14 ks. 



TABLE 4 

Spectral Fits to M15 X-3, 2004 CXO 



Model 


Nh X 10^0 


r 


Xl/dof 


Lx 


Lx,NS 


kT 


Pow 

Pow+NSATMOS 
MEKAL 


4.6* 
4.6* 
4.6* 


1.49±0.13 
1 42+0-18 


0.89/29 
0.90/28 
0.96/29 


5.9^.^6 X 10^^ 
5.8ig-J X 1033 
5.9lg;^ X 1033 


1.7it;? X 10^^ 


09+26 
14ir 



Note. — Power-law, power-law plus neutron star atmosphere, and mekal model fits to spectrum of M15 X-3 in 2004. Errors are 90% 
confidence for a single parameter. X-ray luminosities for 0.5-10 keV, in ergs/s. kT in eV for NSATMOS model, keV for MEKAL model. 
*: fixed to the assumed cluster value. 



